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ABSTRACT 



Aims. We investigate the nature of the source GDS J0332 18.92-275302.7 at redshift ~ 5.56. 

Methods. The spectral energy distribution of the source is well-sampled by 16 bands photometry from UV-optical (HST and VLT), 
near infrared, near infrared (VLT) to mid-infrared (Spitzer). The detection of a signal in the mid-infrared Spitzer/IRAC bands 5.8, 8.0 
fim - where the nebular emission contribution is less effective - suggests that there is a Balmer break, the signature of an underlying 
stellar population formed at earlier epochs. The high-quality VLT/FORS2 spectrum shows a clear Lya- emission line, together with 
semi-forbidden Niv] 1483.3-1486.5 also in emission. These lines imply a young stellar population. In particular, the Niv] 1483.3- 
1486.5 feature (if the source is not hosting an AGN) is a signature of massive and hot stars with an associated nebular emission. 
Conversely, it may be a signature of an AGN. The observed SED and the Lya emission line profile were modeled carefully to 
investigate the internal properties of the source. 

Results. From the SED-fitting with a single and a double stellar population and from the Lya modeling, it turns out that the source 
seems to have an evolved component with a stellar mass of ~5xl0 10 M Q and age ~ 0.4 Gyrs, and a young component with an age of 
~ 0.01 Gyrs and star formation rate in the range of 30-200 M Q yr~ [ . The limits on the effective radius derived from the ACS/z850 and 
VLT/Ks bands indicate that this galaxy is denser than the local ones with similar mass. A relatively high nebular gas column density is 
favored from the Lya line modeling (N HI >10 2[ cm~ 2 ). A vigorous outflow (~ 450 km/s) has been measured from the optical spectrum, 
consistent with the Lya modeling. From ACS observations it turns out that the region emitting Lya photons is spatially compact and 
has a similar effective radius (~ 0.1 kpc physical) estimated at the ~1400A rest-frame wavelength, whose emission is dominated by 
the stellar continuum and/or AGN. The gas is blown out from the central region, but, given the mass of the galaxy, it is uncertain 
whether it will pollute the IGM to large distances. 

We argue that a burst of star formation in a dense gas environment is active (possibly containing hot and massive stars and/or a low 
luminosity AGN), superimposed on an already formed fraction of stellar mass. 

Key words, galaxies: formation — galaxies: evolution — source GDS J033218. 92-275302. 7 



1. Introduction 

In the past few years, dedicated space-borne and ground- 
based observatories and refined techniques have allowed us 
to discover and analyze galaxies at increasingly large dis- 
tances. It is common practice in observational cosmology 
to select efficiently star-forming galaxies (e.g. Lyman break 
galaxies, LBGs, or Lyman alpha emitters, LAE s) and active 
galac tic nuclei (AGN) u p to red s hift 6.5 (e.g. ISteidel et all 
1999; Dickinson eTaR l2004t iGiavalisco et all 12004a: 



Bouwens et al. 



2006; Cristiani et 



20061; Taniguchi et al. 
al.L 12004 iFanetal.L 



20051; lAndoetal. 



2003 
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Fontanot et al. 

2007) or spheroidal/fossil massive galaxies up to redshift 4, 
whose apparent morphology has recently i ntroduced a sub class 
of ul t radense objects (e.g. iDaddi et al. L 12005b ICimatti et al. L 
120081; Ivan Dokkum et alll2008t|Bui trago et al. L I2008I) . 

High-redshift galaxies have been studied through deep multi- 
wavelength surveys, with the aim of maximizing the information 
on the energetic output of the sources spanning a wide range of 
the electromagnetic spectrum from X-ray to radio wavelengths. 
This approach has shown its efficiency of constraining lumi- 
nosity fun ctions up to redshift 6-7 an d down to few percent 
of L* (e.g. lBouwens et al. 1 (120071 12008)). stellar mass functions 
(for m asses higher than \0 9 M a ) up to redshift 6 (Fontana et al. \ 
2006 . 120091; lEvles et al. L 120051 120071: IStark et all 120071 1200' 



lYan et al. L 120051 : lLabbe et al. L 20061) and morphology evolution 
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(iTaniguchi et all l2009t IConselice et ail 120081: iFerguson et all 
2004)". 

Detailed studies of internal the properties of high-redshift 
galaxies, such as information about hot stars, dust, ionized gas 
in HII regions, and the large-scale outflows of neutral and ion- 
ized interstellar mate rial (ISM) are now b e coming feas i ble up 
to redshift 4 -5 (e g Ishaplev et all l2003t lAndo etal.L l2007t 
lOuchi et al. L 120081: IVanzella et alTl2009l) . Useful information 
from the Lya profile modeling of high-redshift galaxies with ra- 
diative transfer codes is producing interesting constraints on the 
dyna mical and physical stat e of the ISM and ionizing so urces 
(e.g. lVerhamme et al.l(l2008l) : ISchaerer & Verhammel(l2008l) ). 

From this point of view, thanks to the combination of depth, 
area, and multivawelength coverag e, the Great Obs e rvatorie s 
Origins Deep Survey p roject (see iDickinson et al. I d2003b); 
iGiavalisco et all (12004b) for a review about this project) is ideal 
for studing galaxies at high redshift and the connection between 
photometr ic, spectroscop i c, and morphological - size proper- 
ties (e.g. IPentericci et all 120071 l2009t IConselice et al. L l2008t 
120061 r 



iRavind ranath et all |2006[ IVanzella et all 120091). a ndlheir rela- 

007h. 



tion with the environment (e.g. lElbaz et al. L 12007). 

Along with enabling a systematic study of normal galaxies, 
multi-frequency surveys over large areas and depth also allow 
us to discover rare objects. In particular, a new class of objects 
showing prominent N iv] 1486 emission have recently been re- 
ported. Such a feature is rarely seen at any redshift. A small 
fraction (1.1%, 1.7<z<4) of the QSO sample extracted from the 
SDSS fifth data release is nitrogen rich, showing N iv] 1486 or 
Nm] 1750 emission lines and Nv] 1240 als o in emission (typi - 
cally stron ger than the rest of the population dJiang et al. 1 12008). 
Similarly, Gli kman et al. I (120071) discuss the discovery of two 
low-luminosity quasars at redshift ~ 4 with Lye and C iv lines, 
moderately broad N iv] 1486 emission, and an absent N v] 1240 
line. In these particular cases, the blinding intensity of the cen- 
tral engine is reduced, allowing study of the properties of the 
host galaxy. Fosbur v" ;t al.lj20 03) report on an HII lensed galaxy 
at redshift 3.357 (the Lynx arc) whose spectrum shows Niv] 
1486, Om] 1661, 1666, Cm] 1907, 1909, as well as the absence 
of the N v] 1240 line. Their modeling of the spectrum favors a 
hot (T~80000 K) blackbody over a n AGN as the ionizing source. 
Alternativelv. lBinette et al. I d2003) sug gest an obsured AGN as a 
photoionizing source of the Lynx arc.lVillar-Mar tin et all {2004) 
propose a population of Wolf-Rayet (WR) stars as the ioniza- 
tion source for the same HII galaxy, with an age below 5 Myr 
that contributes to a fast enrichment of the interstellar medium. 
In this scen ario the stars i nvolve d are much colder than those 
proposed in lFosburv et alJ d2003). 

In the present work we report on the source GDS 
J033218.92-275302.7 at redshift 5.563 located within the 
GOODS southern field, for which extensive information (pho- 
tometry, spectroscopy and morphology) is available. The galaxy 
has b een discovered during the ESO/FORS2 spectroscopic sur- 
vey (IVanzella et all d2006l) ). We focus our attention on this 
source because it shows several unique characteristics. First the 
high S/N spectrum exibits a relatively bright Niv] 1486 fea- 
ture in emission, a unique example among the more than 100 
spectra of high z starburst ga laxies that were obta ined from 
the GOODS/FORS2 campaign dVanzella et all ((2008)). Second, 
while the spectrum shows a bright Lya line indicating a young 
stellar component, the photometry shows a prominent Balmer 
break indicating that there is also an evolved component. Finally, 
the bright IRAC flux suggests that this is a massive galaxy, espe- 
cially interesting if one considers its very high redshift. 



Inde ed the object wa s alrea dy noted by several authors; for 
example, Font anot et al. I (120071) selected it as an AGN candidate 
on the basis of morphological and color considerations, but dis- 
carded it from th e sample because of i ts peculiar galaxy - like 
optical spectrum. Wikli nd et all (|2008) report it among the 11 
candidates with photometric redshifts in the range 4.9<z<6.5, 
dominated by an old stellar population, with ages 0.2-1.0 Gyr 
and havin g very high stellar masses, in the range (0.5 - 5) x 10 11 
M . Also Star k et al.l (120071) re port for this source a st ellar mass 
highr than 10 n M Q . Similarly IPentericci etafl d2009l) note that 
this is one of a handful of bright Lya emitting LBGs at high red- 
shift with an evolved population, indicating that not all Lyman 
alpha emitters are young primeval galaxies. 

As noted rece ntly by ISchaerer & de Barros I (120091) and 
Raiter et al. (2009), strong nebular emission lines may bias the 
result of the SED-fitting of high-redshift galaxies (e.g. [Oiii] 
4959-5007, Ha). Indeed the apparent photometric breaks might 
actually be produced in some cases by the boost of some of the 
lines. Something similar may be happening in this case, there- 
fore it is important to quantify the strength of these lines and 
their influence on global photometry. 

In the present work we perform dedicated SED-fitting allow- 
ing single and multiple stellar populations, and important infor- 
mation is extracted from the Lya- profile modeling. Together with 
the morphological appearance, constraints have been placed on 
the stellar mass density, ages, gas, dust content, and outflows. 

The work is structured as follow. In Section 2 a summary 
of the photometric, spectroscopic, and morphological proper- 
ties is given, and in Sect. 3 the possible scenarios are discussed 
about the nature of the source. Section 4 describes the SED and 
Lya modeling, and in Sect. 5 we discuss the results. Section 6 
concludes the work. In the following the standard cosmology is 
adopted (H =70 km/s/Mpc, Q M =0.3,Q A =0.7). If not specified, 
magnitudes are given on the AB system. 



2. Source GDS J03321 8.92-275302.7 

Source GDS J0332 18.92-275302.7 is located within the south- 
ern field of the Great Observatories Origins Deep Survey. The 
multi-wavelength observations consists of deep U, R (VLT), 
#435 , V 60 6 ' '775 and Z85o (HST), Js, H, Ks (VLT), 3.6, 4.5, 5.8, 
8.0, and 24 ymi (Spitzer) bands. Moreover observations in the X- 
ray and radio domains are available from Chandra and the Ver y 
Large Array, respectively dLuo et all 120081: iMiller et al. 1 12008). 

A considerable part of the spectroscopic information of that 
field has been collected by the VLT/FORS2 spectrograph, which 
has produced about one thousand redshift determinations (with 
a resolution of ~13A, at ~8600A), between redshift 0.5 and 
6.2, in particular more than one hundred LBGs have been con - 
firmed at redshift beyond 3.5 (IVanzella et alll2005ll2006l 12008). 
The VLT/FORS2 spectroscopic survey has been complemented 
in the redshift interval 1.6 < z < 3.5 and at z < 1 by the 
VLT/VIMOS spectroscopic survey, whi ch is producing more 
than 5000 spectrosc opic identifications (Popes so et al. L 2008; 
iBalestra et al. L l2010l) . Such measurements increase the spectro- 
scopic informa t ion available from previous works (e.g. VVDS 
iLe Fevre et all 120051 ISzokolv et al. L l200l . 

Source GDS J0332 18.92-275302.7 with z 850 =24.6 1+0.03 
was selected as a V606-band dropout an d was confirmed to be at 
redshift 5.563 (redshift of the Lyor line. IVanzella et all d2006l) ). 
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Fig. 1. 1 -dimensional spectrum of the galaxy discussed in the 
present work. The Lya and N iv] 1486 emission lines are evi- 
dent, together with the detection of the continuum and the IGM 
attenuation. The inner box shows a zoom of the Lya line, the 
asymmetry and the red tail are visible (the latter is marked with 
a solid segment, 40A long, or 1500 km The transmissions 
of the ACS V606 , '775 . and 2350 are also shown. The exposure 
time was 14400s. 



2.1. UV spectral properties 

The main spectral features are the Lya emission line 
(EW~60A rest frame), the break of the continuum blueward of 
the line, and the semi-forbidden emission line N iv] 1486, a dou- 
blet AA 1483.3-1486.5A (see Fig. Q] and TableQ]for a physical 
quantities summary). The detection of the Niv] 1486 in emis- 
sion is unusual fo r LBGs. However, th is atomic transition has 
been identified by iFosburv et al.l d2003l) i n the Lynx arc and in 
a sub class of OSOs (e.g JGiikman et all (120071) . Baldwin et al. I 
(120031) ). In the following the main properties of the UV spectrum 
are described: 



As shown in the top panel of Fig. [2] an emission line at A 
9742A is detected. The spectral resolution is in principle suf- 
ficient to resolve the double profile of the two Niv] 1486 
components, i.e 1483. 3A and 1486. 5A. We interpret this 
feature as the detection of one of the two components. A 
first possibility is that this line is the 1483. 3A component: 
in this case the redshift turns out to be higher than the ob- 
served Lya redshift. The Lya emission from LBGs is com- 
monly observed to be redshifted relative to the systemic ve- 
locity traced by other, non-resonant emission lines o r stella r 
absorption lines (IShaplev et al.l d2003l).lTapken et all (2001), 
IVerhamme etaTI (120081) . IVanzella et al.l (120091) ). Therefore, 
it would be unusual if, in this object, the Niv] redshift were 
higher than that from Lya . The other possibility is that we 
are detecting the 1486.5A component at redshift 5.553, sug- 



Table 1. Summary of the physical quantities derived from the 
spectral features and morphological analysis. 

Obs. Spect. properties 



z(Lya 1215.7) 
FWHM ( Lya 1215.7) 
EW ( Lya ) 
L( Lya ) 
SFR Q (Lya) 



5.563 

600+100 km s- [ 
59 + _\fk (89 A from phot.*) 
3.8*xl0 43 erg «' 
31 Me/yr 



z(Nrv] 1486) 
FWHM (Niv] 1486) 
EW (Nw] 1486) 
L(Niv] 1486) 



5.553 (1486.5A) 
400+100 km s- 1 
22+^A (33A from phot.*) 



L(Nv 1240-1243) 



1.3_^xl0 43 



erg sec ') 



<1.06±0.57 xl0 4I erg sec' 1 



L(Sirv 1393.8-1402.8) 



<1.13±0.73xl0 42 er£ sec--' 



L(Civ] 1548.2-1550.8) 



<6xl0 42 erg sec-- 1 {2a limit)t 



V(Lya-Niv] 1486) 



+457±50/tm s 



Morph. (rest-frame) 



(kpc physical) 



r e (Lya) (GALFIT) 
r e (1400 A) (GALFIT) 
r e (3300A) 

S/G (1400A) (SExtr.) 
area (1400 A) (SExtr.) 



0.08±0.01 kpc (ACS ;' 77 5 band) 

0.1 1±0.01 kpc (ACS Z850 band) 

< 0.9 kpc (0.6" seeing, ISAAC Ks band) 

0.83 (S/G, ACS zgso band) 

8.8 kpc 2 (AREAF 303 pix., z 850 band) 



f continuum not detected, 
★ adopting the continuum 



2<r of the noise fluctuations is adopted, 
derived from the ACS zgso band. 



gesting a high-density limiQ and a velocity offset (i.e. a pres- 
ence of an outflow) between Lya and Niv] lines of +457 
km/s (dz=0.01). A further proof of this possibility is that the 
redshift and outflow estimated from the spectrum are consis- 
tent with the results from the Lya profile modeling discussed 
below (see Sect. 4.2). Therefore, in the following we assume 
the line to be the 1486. 5A component. 

2. There is no detection for the SilV 1393.8-1402.8 doublet, 
either in emission or absorption, and similarly for the N v 
1240-1243 doublet; their luminosity limits are reported in 
TableQ~](see also Figure|2]i. 

3. From the 2-dimensional spectrum, the FWHM of the spa- 
tial profiles of the Lya and the continuum (by collapsing 
columns along the Y-axis, see Fig. are fully comparable, 
~ 0.7 arcsec. This is also compatible with the seeing dur- 
ing observations, 0.7 arcsec. More interestingly, a better con- 
straint on the Lya extension comes from the ACS (775 band. 
As shown in Fig. Q~| the (775 band is mainly probing the UV 
emission region between the Lya line and 1300 A (less than 
100 A rest-frame), and the part blueward of the Lya line is 
strongly attenuated by the IGM absorption. Since the Lya 
equivalent width is ~ 60A, it turns out that the ACS (775 im- 
age is dominated by the Lya emission. In Sect. 2.3 we show 
that the morphological properties derived from the (775 band 
(probing the Lya line) and the zsso band (not containing the 
Lya line and probing the emission at 1400A) are similar, e.g. 
the effective radii of the two are the same order. This implies 
that the spatial extension of the Lya line is similar to the 
emitting region at 1400A rest-frame. 

Figure [3] shows the contour plots of the 2-dimensional Lya 
region. Various sigmas above the noise fluctuation are reported 
from 1 to 20. It is evident that the asymmetric shape in the wave- 
length domain (that extends at least ~ 40A (3<x)). If the Lya 



1 the ratio of those components is re lated to the electron density, 
iRaiter et alJ d2009h . lKeenan et al. I dl995h 
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Fig. 2. Three zoomed regions of the ID spectrum (red dotted 
line is the solid black line smoothed over two pixels). Dotted 
plots show the rescaled sky spectrum. Filled circles mark the 
position of the lines at the redshift 5.553 (case in which only 
the N iv] 1486.5 component is detected). Dotted open circles de- 
note the other (but less probable) case in which the Niv] 1483.3 
component is detected (z=5.568). Red filled triangles mark the 
positions at the observed Lya redshift, z=5.563. Top: the posi- 
tions of the Niv] 1486 and Civ] features are shown. In the sky 
spectrum two sky emission lines are marked: (1) 10082.46A and 
(2) 10124. 01 A (also shown in the 2-dimensional spectrum of 
Fig. [8}. Middle: the region where N v is expected to lie is shown. 
The open blue triangle indicates the estimated intrinsic Lya red- 
shift from the Lya modeling. Bottom: The region around the 
SilV feature is shown. In the middle (/l~8150A) and bottom 
(/l~9170A) cases, the sky-lines contamination is minimal. 




Wavelength (A) 

Fig. 3. Contour plot of the 2-dimensional Lya region, spa- 
tial direction (pixels, 1 pixel=1.5 kpc) versus wavelength (A , 
20A corresponds to 751 km s^ 1 ). The spectral interval of 7960- 
8280A and spatial extension of 22.5 kpcs are shown. Thin dot- 
ted and continuous black lines mark Ox, 5x, lOx, 20xcr above 
the mean signal. Thick lines mark lcr (red) and 3<x (blue). 
Horizontal dotted lines indicate the FWHM of the profiles (spa- 
tial extension) by collapsing columns in the Lya region and on 
the continuum region. 
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Fig. 4. Cutouts of the source GDS J033218.92-275302.7. From 
left to right: the S435 , y 606 , ins , and z 850 HST/ACS 
bands. The box side is 1.0 arcsec (6 kpc proper at the 
redshift of the source). Drawn from the V2.0 ACS catalog 
http://archive.stsci.edu/prepds/goods/ The last right box is the 
ISAAC/Ks band. The box side is 3 arcsec. 



emission arises from a simple expanding shell of material, then 
it is expected to have less of velocity width at its outer extremes 
than along the line of sight to the central region, where the trans- 
verse velocity component tends to zero. 

2.2. Photometric properties 

Fig. |5] shows the overall SED (black squares) and Table [2] 
summarizes the multi-band photometry of the source collected 
from different instruments mounted on ground and space-based 
telescopes. Magnitudes, errors and 1-cr lower limits (l.m.) 
are derived from the MUSIC catalog, (Grazi an et al. 1 12x1 
ISantini et al.L 120091) . There are other two WFI U bands obser- 
vations not reported in the table, U35 and U38, with a slightly 
different filter shape. Their lower limits are 27.84 and 26.75, re- 
spectively, much shallower than the limit provided by the VLT 
U vimos ■ 

Despite the high redshift of the source discussed here, its 
(2775-Z850) color is 0.59, significantly bluer than the typical 
threshold of ~1.3 adopted to select galaxies beyond red shift ~ 
5.5 (e.g. iDickinson et ail l2004t iBouwens et aU 120071) . This 
comes from the contribution of the Lya emission to the flux in 



the z'775 filter (see Fig. [TJ, that decr eases the (ins-z x w) co lor by 
about 0.7 magnitudes (as shown in IVanzella et al.l d2009l) ). For 
this reason it has been selected as a V606-band dropout source. It 
is also an R-band dropout source if referred to the ground-based 
photometry, see Table [2] As shown in Fig. [4] the source has been 
detected in the Va)6 band with a magnitude of 27.63+0.22, show- 
ing an attenuation of ~ 94% to respect the 1400A emission (z&so 
band), which is c onsistent with th e average IGM transmission at 
this redshift re.g. lSongailall2004l) . 

Apart from the IGM attenuation that influences the bands 
bluer than the zsso, the main feature of the SED is the dis- 
continuity detected between VLT/ISAAC J, H, and Ks bands 
(A < 3400A rest-frame) and the Spitzer/IRAC channels (A > 

5400A rest-frame), see Fig. [5] Such a discontinuity is con- 
sistent with estimates available in t he lite r ature for the sam e 
object (FIREWORKS IWuvts et al.l (l2008h . IStark et alJ (120071) . 
IWiklind et all (l2008h . lRaiTer et al.l (l2009h ~ 

The typical uncertainties (lcr) span the range 0.03 to 0.11 
going from the HST/ACS, VLT/ISAAC and Spitzer channels 1 
and 2, for the last two Spitzer bands (5.6 and 8 pm) the errors 
increase to ~ 0.2/0.3 mags. 
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Fig. 5. The resulted template fitting over the MUSIC multi-band catalog of the GDS J033218.92-275302.7. Left: Single stellar 
population modeling. Blue solid line is the fit adopting the maximum ratio [O ni]/[0 n] (prescription "Single/[0 in] max" in TableO, 
and red dashed line is the fit with the Schaerer & de Barros (2009) method ("Single SB09" in the same table), see text for details. 
Right: Double stellar population modeling. Green dashed line shows the evolved component (age ~QA Gyr), the thin cyan line the 
young contribution (age ~ 0.01 Gyr), and black thick line the best fit composition of the two. 



Table 2. Summary of the photometric information (magnitudes and 1-cr errors) for our source. 



U VIMOS t 


-6435 


^606 


R 


'775 


ZS50 


Js 


H 


Ks 


3.6//ra 


4.5/j.m 


5.8//m 


8.0yt/m 


24fim 


29.49 


29.86 


27.86 


27.11 


25.22 


24.61 


24.68 


24.80 


24.21 


22.82 


23.00 


23.37 


23.18 


22.61 


lcr l.m 


lcr l.m 


0.179 


0.195 


0.039 


0.031 


0.116 


0.151 


0.140 


0.014 


0.024 


0.235 


0.198 


lcr l.m. 



t Other two WFI U bands are available, U35 and U38 with slightly different filter shape. The lower limits are 27.84 and 26.75, respectively, 
much shallower than the limit provided by the Uvimos ■ 



2.3. Morphological properties 

Image cutouts of the isolated source GDS J033218.92-275302.7 
in the B435 , V 606 , '775 , Z850 (HST/ACS), and Ks (VLT/ISAAC) 
bands are shown in Fig. [4] where each box is 1.0 arcsec wide 
for the ACS figures (6 kpc proper at the redshift of the source), 
while it is 3 arcsec on a side for the ISAAC Ks image. We con- 
sidered the ACS/Z850 and the ISAAC/Ks bands to derive basic 
morphological quantities (see Table[TJ. 
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The uncorrected-PSF effective radius (Re) available from the 
ACS/GOODS public catalog v2.0 is 3.37 pixels. The same 
quantity derived from a sample of 80 stellar-like sources 
(SExtractor star/galaxy index larger than 0.97, 1.0=star, 
0.0=extended source) with zm) magnitude in the range 23.5- 
25 gives a median value and lcr percentiles of 2.590^'^, 
implying that this source is not a stellar-like object (see 
Fig. |3J, The same result is obtained for the 2775 band (~6<x 
from the median value of the stars). We note that the 
SExtractor star/galaxy index of the source is quite high, 0.83, 
but lower than the typical value of the stars. Therefore, even 
though the present galaxy is clearly a compact source in the 
UV rest-frame, it is marginally resolved both in the 2775 and 
z 850 bands. 



Fig. 6. The effective radius measured with the SExtractor algo- 
rithm (pixel units) as a function of the z&so magnitude for a sam- 
ple of stellar-like sources selected having S/G classifier larger 
than 0.97 (crosses). Solid line is the median value and the dotted 
lines the lcr percentiles. The open red circles are stars confirmed 
spectroscopically. The filled blue circle indicates the source de- 
scribed in the present work. 



To derive PSF-correcte d morphological pa rameters, we ran 
the GALFIT program dPeng et al. 1(120021) ) in both bands. 
The morphological shape of the source is not particularly 
complicated so a good fit is reached (reduced ;f 2 = 0.591) 
by adopting a simple Gaussian profile (Sersic model with 
n=0.5) and leaving the Re, the axis ratio B/A, the coordi- 
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Fig. 7. The ACS zg 50 image of the source (left, 100x100 pix- 
els) and the residual image after subtraction of the model de- 
rived from GALFIT (middle). In the right panel the behavior of 
the residuals (standard deviation calculated on a 20x20 pixels 
area centered on the source) as a function of the Re (fixed dur- 
ing GALFIT runs) is shown. The minimum value corresponds 
to Re ~0.65 pixels, i.e. 0.02 arcsec or 0.12 kpc proper. The dot- 
ted line indicates the median value of the background residuals 
(standard deviation calculated in the blank regions). 



nates X, Y, the magnitude, and the position angle as free pa- 
rameters (SExtractor estimates were used as a first guesses). 
In the left panel of Fig. [7] the zxso band image of the galaxy 
is shown, and the residuals map provided by GALFIT, as a 
result of the subtraction of the best-fit model from the origi- 
nal galaxy, do not show significant structures (middle panel 
of the same figure). An effective radius /?e=0.62±0.05 pix- 
els (0.11+0.01 kpc physical) and B/A=0.61+0.14have been 
obtained. We explored how the variation in Re affects the 
residuals fixing n=0.5 and B/A=0.61. GALFIT was run 80 
times by varying the Re from 0.025 to 2.0 pixels with a step 
of 0.025. Its behavior is shown in the right panel of Fig. [7] 
where the minimum of the residuals for Re ~0.650 is clear, 
still with acceptable values between Re 0.3-1.0 (pixels). It is 
worth noting that, if we perform the same fit in the (775 band, 
the best Re is slightly smaller than the Z&50 band (0.425 pix- 
els, i.e. 0.013 arcsec or 0.08+0.01 kpc physical), implying 
an even more compact region for the Lya emission. 
- Ks band (A rest ~ 3300A). 

To approach the optical rest-frame wavelengths, we exploit 
the information derived from the near infrared (NIR here- 
after) observations (VLT/ISAAC Ks band). The resolution 
from the ground is not comparable to that obtained from 
the space; however, we note that the present object has an 
FWHM fully consistent with the seeing at the epoch of 
observations of that particular region of the GOODS field 
(0.50 arcsec), i.e., it is not resolved. Assuming a Gaussian 
shape with an FWHM of 0.50 arcsec, it turns out that the 
r<,[3300A] is 0.9 kpc physical. For the Gaussian profile the 
i?e=(FWHM/2)/1.738, i.e. the radius containing half of the 
total light. Since the source is not resolved, the derived Re is 
an upper limit at this wavelength. 



There are three other sources (with secure redshift, qual- 
ity 'A") in the FORS2 sample that show a similar z&so UV 
morphology of the source discussed here (with a SExtractor 
Re smaller than 3.5 pixels (0.105 arcsec)) . One is the QSO 
GDS J033229.29-275619.5 (iFontanot et al. I (I2007h 1 fully com- 
patible the stellar PSF and two other LBGs: GD S J033217.22- 
274754.4 at redshift 3.652 already described in Vanzell a~et alJ 
(2008) with a double structure of the Lya profile and 
GDS J033240.38-274431.0 at redshift 4.120, which shows a 
continuum with clear Lya emission and no other features. 
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Fig. 8. Extracted VLT/FORS2 2-dimensional spectrum of 
the galaxy discussed in the present work (the spectral inter- 
val 9403-10164A is shown). As a check of the wavelength 
calibration the skyline position derived from the FORS2 spec- 
trum are reported (the sky line measurement performed by 
VLT/UVES and Keck/HIRES are indicated within parenthesis 
www.eso.org/observing/dfo/quality/UVES/pipeline/skyjpectrum.html). 
The expected positions of the C iv] feature are also shown with 
thick arrows: case (A) at the redshift of the observed Niv] 1486 
(z=5.553) and case (B) at the redshift of the Lye line (z=5.563); 
see text for details. 



3. Possible scenarios for GDS J03321 8.92-275302.7 

3.1. A chance superposition ? 

If we interpret the emission line detected at /l=9742A as a [0 11] 
3727 foreground emitter superimposed on a background and 
brighter higher-z (z=5.563) LBG, its redshift would be 1.614. 
It is well known that in the GOODS-S fiel d there is an overden- 
sity structure at that redshift, z=1.61 (e.g . IVanzella et all [2008; 
ICaste llano et al. L 120071: iKurk et al.U2009l) . However we can ex- 
clude this possibility on the basis of the ACS morphology, shows 
a compact and circular shape (see Fig. |4j and from the ultradeep 
U band observations carried out by the VLT/VIMOS instrument 
dNonino et alJ d2009l) ). provide an upper limit of ~30 AB at lcr 
(also it has not been detected in the ACS B435 band image). The 
source should be detectable in the blue if there is star formation 
activity traced by the [On] 3727 line. Moreover, assuming a flat 
continuum at that lcr limit (30 AB), the rest-frame [On] 3727 
equivalent width would be larger than 10 4 A. Even though we 
consider this possibility largely unlikely, we note tha t an exam- 
ple of strong [O n] 3727 emitter has been reported bv lStern et al. I 
d2000b . 

3.2. Is it an AGN ? 

The spectral range up to 10164 A allows us to detect possi- 
ble emission lines testing for the presence of an AGN (e.g. 
Nv] 1240-1243, SilV 1393.8-1402.8, and Civ] 1548.2-1550.8). 
Those features are routinely detected in spe ctra of the most ob- 
scured AGNs fe.g. lPolletta etIHl2006ll2008h . 

3.2.1. Line emission in the UV 

- Nv 1240-1243 feature. 

Observationally, the Nv emission is often present in the 
AGN case. The FORS2 spectrum allow us to measure emis- 
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sion features down to ~2xl0 erg/s. No N v line is detected 
(see Fig . |2j - 
- Civ 1548.2-1550.8 feature. 

Unfortunately, the Civ line is at the very red limit of the 
observed spectrum, which in this particular slit position 
is determined by the detector end. Figures [8] and [2] (top 
panel) show the 2-dimensional and the 1 -dimensional spec- 
tra zoomed at the red edge, respectively. On one hand, 
if the redshift of the Civ feature is higher than 5.565 
(z(Civ)>5.565, higher than the Lya redshift), then the dou- 
blet is completely out. Alternatively, if z(C iv)<5.554, then 
the doublet is completely in. In general, LBGs show the ob- 
served Lya redshift as higher than the other spectral fea- 
tures because of its asymmetry which typica lly arises from 
backscattering of the reced ing material (e.g. IShaplev et al.L 
l2003t IVanzella et ail l2Q09h . Similarly, in the AGN case, a 
blueshift of the Civ feature with respect to the observed 
Lya peak is typically observed, of several hundred kilo- 
meters per second. In particular a blueshift of ~ 600±100 
km with respect to the Lya line is measured in the SDSS 
QSO composite spectrum (e.g. see Table 4 and Figure 9 of 
IVanden Berk et al. I d200lh l. If this is the case, the C iv fea- 
ture should fall in the available spectrum, and its luminosity 
limit is ~6xl0 42 erg/ s at 2cr of the background noise fluctu- 
ation. 

The limits and luminosities estimates on Lya , N v, Si iv, and 
C iv are reported in Table Q] 

3.2.2. X-ray emission 

We further note that this source has not been detected in the X- 
ray by the 2 Ms Chandra ultra-deep observations, neither in the 
MIPS 24/mi by Spitzer (with lcr lower li mit of 22.61 AB) n or 
by the VLA at 20 cm down to 8 fiJy (at lp-. lTozzi etaTI (120091) 1 . 

A correlation between 2-10 keV X-ray luminosity and [O m] 
5007 or Ha e mission line luminosit ie s is observed for galaxi es 
at z < 1 (e.g., iPanessa et ail (120061) . ISilverman etaD (120081) 1 
Assuming that this correlation also holds at higher redshifts, then 
we may use the X-ray luminosity to derive a constraint on the 
fluxes for AGN-powered emission lines in the IRAC bands. 

From the current 2 Ms observations, there is no detection 
at the position of t he source ( 3<x limit of ~ 3xl0 43 erg/s at 3- 
13 keV rest-frame. iLuo et~a 1. (2008)). This limit roughly cor- 
responds to an upper limit for both Ha and [O m] 5007 lumi- 
nosities of ~ io 40 - 5 ~ 42 - 5 and 10 41 ~ 43 erg s _1 , respectively. Such 
values are affected by large uncertainties in the assumed rela- 
tions (intrinsic scatter) and the limit derived from the 2 Ms im- 
age. However, if compared to the (at least) one magnitude jump 
between the VLT/ISAAC and Spitzer/IRAC magnitudes, these 
estimations suggest that, besides a line contribution to the IRAC 
magnitudes, there is also a significant contribution from stellar 
emission beyond 5000A rest-frame, i.e., of a relatively evolved 
stellar population (see Sect. 4). For example a line luminosity of 
7xl0 45 erg s _1 is needed to boost the 4.5fim AB magnitude from 
24 to 23 (adopting a bandwidth of lOlOOA jFazio et al. I d2004 ). 

3.2.3. A rare class of QSOs: an open possibility 

It is interesting to compare our N iv] 1486 emitter spectrum with 
the composite spectra of quasars available in the literature. This 
has al r eady been done by se veral authors (e.g. [Baldwin et al. I 
d2003l) . lGlikman et al. I d2007l) . and lJiang et alTd2008l) ). None of 
the published average quasar spectral templates show any trace 



of N iv] 1486 emissi on. Neverthel e ss, f ocusing the attention on 
this spectral feature, iBentz et aI71 d2004l) compile a sample of 
6650 quasars in the range 1.6<z<4.1 showing the Niv] 1486 
line (othe r than the N ml 1 750), and more recently an updated 
work by Jiang et al. I d2008) (on SDSS data release 5) reported 
that such objects are ~ 1.1% of the total SDSS quasar sample. 
They also note that for this small fraction, the N v 1240 and Lya 
are much stronger than the rest of the population. We recall that 
our source does not show the N v 1240 line. 

More inte r esting ly and similar to our findings, 
iGlikman et al. I d2007l) have discovered two low luminosity 
QSOs at redshift ~ 4 showing Lya , Niv] 1486 and Civ 
1548-1550 emissions, but no detection of N v 1240. In one case 
the equivalent width of the N iv] 1486 is larger than the C iv one 
(240A vs. 91A), while it is the opposite for the other (24A vs. 
91A). Our source has a luminosity of Mi45=-22.1 (AB) and 
shows a clear Niv] 1486 emission with an equivalent width 
of ~ 22A and FWHM ~ 400 km sT 1 . The Lya line shows a 
narrow componen t with a measur e d FW HM of ~ 600 km 
As performed in Glikman et al. ] (12001 . since the blue side 
of the line profile is absorbed, we forced the symmetry in the 
line by mirroring the red side of the line profile over the peak 
wavelength and computed the Gaussian fit. The narrow-line 
component increases to ~ 750 km The broad-line feature 
(indicated with a segment in the innner box of Fig. Q]) gives an 
FWHM of ~ 3500 km s -1 . This would put the source in the QSO 
regime (velocity width larger than 1000 km s^ 1 ). Therefore, 
the present source may be consistent with the interpretation of 
a low-luminosity q uasar in which the ho st starburst galaxy is 
visible (similarly to lGlikman et al. I d2007l )). 

The study of stellar populations of low-luminosity AGNs 
(e.g. low-luminosity Seyfert galaxies, low-ionization nuclear 
emission line regions, LINERs, and transition-type objects, 
TPs) has been addre s sed f or the local Universe (e.g., 
iGonz alez Delg ado et al. I d2004l) ). but this is still a poorly ex- 
plored regime at higher redshift. While it is beyond the scope 
of the present work to explore the link between the coevolution 
of (circumnuclear) starburst activity and the central black hole 
accretion, we simply note that both AGN and star-formation re- 
quired gas to fuel them, and it happens on different temporal and 
spatial scales, on sub-parsec and typically above few hundred 
parsecs (up t o sever a l kilo-parsecs) regi ons, respectively (e.g., 
iDavies et al~1 d2007l) . IChen et al. 1 d2009l) ). In the present case, 
the size of the UV emitting region is compact, but still resolved 
in the z&so ACS image (as shown in Sect. 2). In summary, the 
presence of an AGN - in a rare evolutionary stage - may be in- 
dicated by the Niv] 1486 and broad Lya features, even though 
Nv 1240, Si iv 1394-1493, and (possibly) Civ 1548-1550 are 
not detected. 

3.3. A multi-burst galaxy in a peculiar stage of evolution ? 

The so urce GPS J033218.9 2-275302.7 has already been ana- 
lyzed in Wikli ndet al.l (12008). who classify it as a "pure" balmer 
break galaxy (their ID #5197). The discontinuity detected be- 
tween the Ks and 3.6fim bands is interpreted as a signature of the 
Balmer break, suggesting a relatively evolved age of stellar pop- 
ulations with a significant stellar mass already in place (age of 
~0.7Gyr an d M*~7-8xl0 10 M p)). A similar conclusion has been 
reached by IStark et al.1 (120071) . who find an even higher stellar 
mass of 10 n M o (their ID 32.8020). 

However, most probably the observed (Ks-3.6fim) color is 
contaminated by emission lines in the 3.6fim band, e.g. [Oin] 
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4959-5007. A similar boost to the flux in the 4.5/^m band may 
come from the Ha emission line. It was also selected as Ha 
emitter by R. Chary et al., private communication. Apart from 
the evident Lya emission, which implies the presence of young 
(< 10 Myr) stars - i.e., some "current"/ongoing star formation 
- significant nebular emission is also robustly supported by the 
detection of the Niv] 1486 line. As mentioned above, a simi- 
lar feature has been identified in the Lynx arc and may indicate 
a short powe rful starburst in whic h very hot and massive stars 
(T>80000 K, i Fosburv et all (120031)") or cooler Wolf-Rayet stars 



are involved ( Villar-Martin et al. 



2004)). A similar blackbody 



ionizing source may be present in this source. The ongoing star 
formation activity would also be responsible for the measured 
outflow, whose spectral signature is in the red tail of the Lya pro- 
file (see Lya modeling in Sect. 4.2). It is beyond the scope of the 
present work to model the ionizing source; nevertheless, we note 
that in a "pure" nebular scenario, the continuum is practically 
flat, and the observ ed "breaks" a re pro duced by strong nebular 
emission lines (see iRaiter etafl (120091) for a dedicated discus- 
sion). Alternatively, a different interpretation suggests a contri- 
bution from a relatively evolved stellar population that produces 
the Balmer break signature (see next section). Given the current 
spectroscopic and photometric information, the following mixed 
scenario may be possible: 1) ongoing active star formation in an 
////-like region that produces nebular emission, as probed by the 
Lya and Nrv] 1486 features, and 2) an already evolved popula- 
tion of stars formed at higher redshift, as probed by the signal 
detected in the IRAC bands, in particular, redwards of the 4.5^/m 
(beyond ~7000A rest-frame). 



4. SED and Lya modeling 

We cannot definitively distinguish between the two scenarios de- 
scribed above, in particular for the explanation of the N iv] 1486 
feature. In either case, even though the source reflects an early 
stage of coevolution of the (circumnuclear starburst) galaxy with 
its AGN or it is an HII source, the features of the host galaxy are 
detected and can be investigated. We therefore need to model the 
SED allowing for multiple stellar populations. Moreover, valu- 
able information can be derived from the Lya line modeling. 
This is performed in the following sections. 

4.1. Modeling the SED 

The SED modeling was performed adopting t he multiwave- 
length GOODS-MUSIC photometric catalog dGrazian et alJ 
(200 6])), and the spectral f itting technique was developed in 
Font ana et aL_ d2003l |2006|) (simil a r to those adopted in other 
works. e.g.. lDickinson et al. I d2003h : lDrorv et al. I d2004l) ). In the 
previous section we pointed out that this galaxy likely contains 
a mixed stellar population, both young stars, as implied by a 
bright Lya line, the Niv] 1486 line, and old stars that produce 
the 4000A break, clearly observed in broad band photometry. 
We therefore model the SED of the galaxy both with a single 
stellar population and with a more plausible double stellar pop- 
ulation. More complicated mixes of multiple stellar populations 
becomes unconstrained by the data given the many degrees of 
freedom of each population. We actually reduce the number of 
degree of freedom by imposing the requirement that both popu- 
lations (old and young) are affected by the same dust extinction 
(with a Calzetti or a Small Magellanic Cloud extinction curve). 
Although this might not be true, we feel that it is plausible for 
very compact objects such as the one we are studying. We fix the 



rest-frame equivalent width of the Lya line to be EW=60A, as 
measured from the spectrum and regardless of the star formation 
rate. 

As discussed in the previous section, the contribution of neb- 
ular lines to the photometry may affect the IRAC magnitudes of 
the 3.6/um and 4.5/vm bands. In particular, the [Om] emission 
contributes to the 3.6 pmn channel and the Ha line to the 4.5 /urn 
channel. While the Ha line can be modeled relatively easily and 
its luminosity can be assumed to be proportional to the global 
SFR through the well-known K ennicutt relations , dise ntangling 
the [Om] contribution is harder. iMoustakas et alJ d2006l) investi- 
gated the [O m] nebular emission line as a quantitative SFR di- 
agnostic and conclude that the large dispersion in the [O m]/Ha 
ratio among star-forming galaxies precluded its suitability for 
SFR studies. 

We therefore treat the [Om] contribution in three different 
ways: 1) by assuming a mean [Om] flux as inferred in local star- 
burst galaxies (corresponding to a ratio f([Oin])/f([On])=0.32), 
2) by assuming a maximum [Om] flux in the 3.6 ym band 10 
times larger than in the previous case, corresponding to the max- 
imum observed [O m]/Ha in star-burst galaxies, and 3) by ne- 
glecting the 3.6 yum band in the fit. 

In the SED-fitting computation, the formal errors of the ob- 
served magnitudes have a minimum value permitted for each 
band. This was done to avoid over- fitting in the^ 2 minimization 
procedure, and it affects only the 3.6 and 4.5 [im bands, whose 
errors are increased to 0. 1 (the minimum permitted) during the 
fit. 

The results of the various fits are reported in Table[3j for each 
model, we report the best-fit (bf) total stellar mass, age, t (the 
star formation e-folding timescale), current SFR, and E(B - V) 
(indicated with EBV) for the single and double populations, as 
well as the minimum and maximum values allowed by the fit (at 
lcr). For the double stellar populations the best-fit ages of the 
young and the evolved components are reported (in the last two 
columns). 

4.1 .1 . A single stellar population model 

Although the single stellar population model with a declining 
exponential SFR is clearly an oversimplification it can set useful 
limits. From Table [3] we see that in all cases the best-fit stellar 
mass is well above 10 10 M© and the age more than 700 Myrs, 
implying a formation redshift z>13. The variation in the [Om] 
flux of a factor 10 does not have a strong impact on the values of 
mass, age and SFR. Even neglecting the 3.6//m band, the stellar 
mass is set to 4.5xlO lo M (with a minimum value of 3.4), but 
the most notable change is that some dust extinction is allowed, 
with a best fit E(B - V) stars =0.06. 

It is worth noting that the stellar mass and ages we find 
are comparable to those derived by Wikl ind et al.1 d2008l) and 
I Stark et all d2007h if the nebular line treatment is not inserted. 
If we insert the [O m] prescription, our estimates become some- 
what smaller, even though still significatively large given the 
redshift of the source (co rresponding to ~ .91 Gyrs after the 
Big-Bang). In particular, IStark et all d2007l) derive a mass of 
1.4 x 1O H M but without including the 5.8yum and 8//m IRAC 
bands in the SED fitting, while Wiklind et al. report a mass of 
7xlO 1() M but assume a photometric redshift of 5.2. 

For comparison we have also fitted the SED with the method 
of Schaerer & de Barros (2009) allowing for numerous emis- 
sion lines. The results obtained (see Table [3] and Fig. [5]) are 
compatible with the two other approaches used here (see also 
next section). In addition, these models predict an intrinsic Lya 
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Table 3. SED modeling: results. 



MODEL 




M bf 




S FR m i n 


SFR bf 


S FR max 


EBV mi „ 


EBV bf 


EBV max 


T 


r 




10 10 M e 


10'°M G 


10 I0 M o 


M Q yr- 1 


M e yr 1 


M Q yr- 1 








Gyrs 


Gyrs 


Single / [O m] ave 


5.2 


6.8 


9.1 


22 


29 


30 


0.00 


0.00 


0.00 


0.71 


0.3 


Single / [O m] max 


5.2 


6.8 


9.1 


22 


29 


30 


0.00 


0.00 


0.00 


0.71 


0.3 


Single / no 3.6 


3.4 


4.5 


6.8 


22 


34 


50 


0.00 


0.03 


0.06 


0.79 


0.6 


Single SB09 t 


4.1 


6.4 


8.4 


17 


36 


41 


0.00 


0.05 


0.10 


<1 


0.7 


Double / [0 m] ave 


4.0 


5.7 


6.9 


22 


54 


249 


0.00 


0.03 


0.15 


0.01,0.4 


0.6,0.1 


Double / [0 m] max 


4.2 


6.2 


6.9 


21 


79 


145 


0.00 


0.03 


0.10 


0.01,0.4 


0.6,0.1 


Double / no 3.6 


3.6 


5.5 


6.0 


23 


81 


146 


0.00 


0.06 


0.15 


0.01,0.4 


2.0,0.1 



t SED-fitting according to the Schaerer & de Barros (2009) method. 



equivalent width of EW(Lya') ~ 49-70 A (1 <x interval), in good 
agreement with the observations. It is worth mentioning that the 
present source is different from those analyzed by Schaerer & 
de Barros (2009), extracted fr om a sample of z~6 star-forming 
galaxies of Ev les et al71 d2007l) . First, the photometric break be- 
tween the NIR bands and the first two Spitzer/IRAC channels 3.6 
and 4.5 pm is at least twice for our source, 1.5 mags. The source 
is brighter in absolute scale, allowing smaller photometric errors 
see Table [2] The main nebular contributors at z>5.6 to the IRAC 
channels 3.6 and 4.5 pm are the H/3, [O m], and Ha lines, respec- 
tively. In this case the H/3 line fall s on th e blue edge of the 3.6 
pm filter transmission (Fazio_etalJ(|2004)), so its contribution is 
further attenuated. However, more importantly, in our case the 
source has also been detected where the contribution of nebular 
lines is less effective (at this redshift), i.e. in the 5.8 and 8.0 pm 
bands. Even though the photometric errors slightly increase in 
these bands, the break is still large ( ~ 1 magnitude). For compar- 
ison, in the sample of lEvles et al~l d2007l) . the source SBM03#1 
(with zg5o = 25.35) has not been detected in the reddest IRAC 
channels 3 and 4. This is the main reason for this galaxy main- 
taining a relatively large estimation of the age and stellar mass, 
even considering the nebular contribution. 



4.1 .2. A mixed stellar population model 

In this case we allow for two stellar populations with different 
redshifts of formation and different star formation histories. The 
only constraint we impose (a part the spectroscopic redshift) is 
that both stellar populations are affected by the same dust ex- 
tinction. In Table [3] the best-fit values for total stellar mass and 
the star formation rate corresponding to the sum of the two con- 
tributions are shown, each taken with the relevant normalization 
factor, while we separately report two best-fit ages for each pop- 
ulation. In all cases, the best-fit solution is made by an older 
stellar component (of ~ 0.4 Gyrs) that contributes to most of the 
stellar mass, while the young component has high values of SFR. 

The global stellar mass remains similar to the previous sin- 
gle population case, regardless of the [O m] treatment, with mean 
values of 5-6 xlO lo M and in all cases greater than 4xl0 10 M Q . 
The SFRs increase by a factor of at least two, and solutions with 
values as high as SFR=250M G yr~ 1 are acceptable. In all cases, 
the e-folding timescale of the star formation rate of the young 
population is close to 100 Myrs. The derived ages of the two 
stellar populations are 0.01 and 0.4 Gyrs, independent of the 
[O in] treatment. The age of the old component is still quite large 
(but less than the previous case), implying a formation redshift 
around z~8. Most important, solutions with small but non negli- 
gible dust extinction are always preferred. 

Formally the best-fit solution among those including all 
bands is the double population with maximum [O ra] contribu- 



tion (x 1 = 4.1). In summary we consider that 5xl0 10 M Q is 
a fair estimate of the total stellar mass and that, considering all 
uncertainties in the data and in the modeling, a solid lower limit 
can be set at ~ 3 x 10 10 M o (the contribution of the young compo- 
nent to the total stellar mass is negligible, ~ 1% in all cases). The 
galaxy contains a stellar population that is at least 400 Myrs old, 
and the average extinction (E(B - y) JM „=0.03-0.06) is smaller 
but not incompatible with the extinction factor that comes from 
the Lya profile modeling (see below). 

4.2. Modeling of the Lya line 
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Fig. 9. The observed Lya- line profile together with the best- 
fitting synthetic profiles described in the text. 



The galaxy shows an Lya emission line with an FWHM of 
600 km s~ l , an evident asymmetric profile, a clear sharp decline 
in flux on the blue side and a red tail of Lya photons extending 
up to ~40A (1500 km s^ 1 ) from the peak of the line (see inner 
box of Fig. [T]and Fig. [3j. 



2 We performed separate fits with enlarged errors in the 5.8 and 
Spin channels to explore their possibile underestimation, and the conse- 
quence is even higher stellar masses, since more weight is given to the 
3.6 and 4.5pm bands. 
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Table 4. Best-fit parameters from the Lya line fitting. Values marked in bold face have been fixed during the fitting procedure. 



Param. 


zsys_V4Q8 


zsys_V5Q8 


zsys 


V4Q0 


V5O0 


free 


n(EI) [cmr 2 ] 


20.8 


20.8 


20.8 


20.8 


21.4 




Vexp [km s" 1 ] 


400 


500 


500 


400 


500 


500!» 01 




3.0 


3.0 


3.0 


3.0 


3.0 


3.0i? 


b [kms- 1 ] 


10.0 


40.0 


40.0 


20.0 


10.0 


10 0+ 150 


z 


5.553 


5.553 


5.553 


5.551 


5.540 


5 540 +!i - 007 


xl 


3.97 


3.22 


3.22 


3.19 


1.785 


1.785 



Lya is a resonance line that undergoes a complicated ra- 
diation transport, with the line formation under the influence 
of numerous parameters: not only dust but also the geome- 
try, k inematics, a nd temperature structure of the neutral ISM 
(e.g. lAhn et al. I (l2003h . IVerhamme etaD (120061) ). These pa- 
rameters influence the line profile and, if sufficient care is 
taken, the line profile itsel f can be used to provi de indepen- 
dent and unique constraints (IVerhamme et al.l(l2008l) ). Using the 
Monte Carlo Lyman-al pha (MCLya) radiation transfer code of 
IVerh amme et al. (2006), we computed a wide array of possible 
emergent line profiles. 

Parameter fitting is performed using a standard least squares 
fitting engine to minimize the^ 2 statistic. Details of the software 
and fitting can be found in Hayes et al. (2009, in prep). The pa- 
rameter space is not entirely unconstrained; e.g., it is possible to 
observationally constrain two of the parameters: the Niv] 1486 
line puts the systemic redshift at 5.553; and the velocity shift be- 
tween Niv] 1486 and Lya constrain the expanding velocity of 
the gas shell to respect the stellar component, V e xp=457 km s _1 
(See Table Q}. Since our grid of shell parameters is discrete, 
we adopted the nearest values of the outflow velocity of 400 
and 500 km s _1 . We ran six independent fits in total for all 
combinations of constraints denoted as follows: zsys_V400 and 
zsys_V5QQ constraining both z and V exp (400 and 500 km/s); 
zsys constraining z; V4Q0 and V500 constraining V exp (400 and 
500 km/s, respectively); and free in which all parameters are fit 
without constraint. 

The results of the fits are presented in Fig. [9] and Table|4] In 
general, all the fits agree with the case free, in particular V500 
and free produce the same values. The Lya modeling favors 
high HI column densities (Nhi > 10 20 ' 8 cm -2 ), outflow veloci- 
ties of 400-500 km/s (consistent with observations when V exp is 
allowed to vary), and a Ti2i6~3.0 for all models with lcr error of 
+ 1/-2 (in the free case), which provides a rough estimate of the 
extinction E(B - V)~0.3+|} 2 . 

We also explored the possibility that the emerging Lya shape 
is caused by a static gas (Ve X n=0). In this case the expected 
double-peaked structure (e.g. IVerhamme et alj d2006l) ) would 
mimic the single peak observed, since the bluer one could be 
self-absorbed (by the galaxy and IGM). From the modeling it 
turns out that high values of the DoppleiQ parameter b and in- 
trinsic FWHM are favored, ~ 160km/s and 1000 km/s, respec- 
tively. This is not surprising since it is a way to drive photons 
away from line center in the absence of a wind. However, in 
all cases the resulting fit worsens in general and, in particular, 
the extended red and wavy tail of the line is no longer repro- 
duced. Conversely, this feature favor the above interpretation of 
backscattered photons from an expanding shell (the presence of 
a wind would agree also with the ongoing star formation activ- 
ity). 



3 Where the Doppler parameter is b = (Vl + V? llrh )°~ 5 as the contribu- 
tion of thermal and turbulent motions. 



It is worth mentioning that we do not expect that the extinc- 
tion undergone by the nebular lines (E(B - V) gas from Lya fit- 
ting) should match what is undergone by the stars (E(B - V) stars , 
from SED fitting). Indeed Verhamme et al. (2008) find a large 
scatter in the relation between the extinction determined from 
the Lya profile fits versus other methods including photometric 
fit and/or measured spectral slopes (see Fig. 12 of their work). 
Calzetti et al. (2001) find an empirical relation of E(B - V) stars = 
0.44 E(B - V)cas, which would make the two results even more 
consistent (within their uncertainties). 

5. Discussion 

5. 1 . Summary of the modeling 

The SED fitting and the Lya modeling indicate that: 

1. The column density of the nebular neutral gas is high, 
A^f/>10 208 . We note that this value is comparable to 
those found for th e damped Lyman-alpha systems, e.g., 
IWoifeetal.l(l2005h . 

2. The outflow velocity derived from the Lya modeling is con- 
sistent with the observed one, and it is relatively high (greater 
than 400 km/s). 

3. A young and an evolved stellar population are both present. 
The first with an SFR in the range 30-200 M Q yr~ l and negli- 
gible contribution to the total stellar mass (1%). The second 
with a stellar mass of ~ 5xlO 1() M and an age of 0.4 Gyr. 

4. The extinctions derived from the different methods are com- 
patible within the 1-cr uncertainties and in general are rela- 
tively small (in the range 0<E(B - V)<0.3). 

Although the signal-to-noise ratios are low, the galaxy is 
strongly detected in both the IRAC 5.8 and 8.0 ymi bands, where 
no nebular lines or strong nebular continuum should contribute 
at z — 5.56. Even if emission lines contribute flux to the shorter 
wavelength IRAC channels, the longer-wavelength IRAC chan- 
nels indicate a significant increase in flux density relative to the 
rest-frame UV continuum, suggesting the presence of a Balmer 
break from an evolved stellar population. 

Summarizing, this galaxy shows several interesting observed 
properties: 1) its stellar mass is still relatively high (M* ~ 
5xlO 1() M ) with a component of already evolved stellar popu- 
lations (-0.4 Gyr); 2) it contains a star-forming component able 
to produce nebular emission lines and with an age of ~ 10 Myr; 
3) a substantial wind is measured both from the optical spec- 
trum and from the Lya modeling, of 450/500 km/s; and 4) the 
source is compact in the rest-frame UV and U-band rest-frame 
wavelengths. 

5.2. An already dense galaxy ? 

As described in the previous sections, the SED fitting analy- 
sis implies a stellar mass of M*~5xlO 1() M , with a significant, 
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evolved component with an age of ~0.4 Gyr. If the very compact 
size measured in the ACS and ISAAC images for the rest-frame 
ultraviolet light can be assumed to apply to the overall distribu- 
tion of the evolved stellar population, then it implies a very high 
stellar mass density: 

1. If we assume a constant size over all wavelengths from the 
1400A to optical bands rest-frame (ACS zsso band, 0. 1 1 kpc 
physical), the stellar density in a spherical symmetric shape 
turns out to be p*=(0.5M*)/(4/3;rr 3 ) ~ 3.5xl0 12 M e kp C -\ 

2. Similarly, assuming a constant size over all wavelengths 
from the 3300A rest frame (ISAAC Ks band, 0.9 kpc physi- 
cal), the stellar density in a spherical symmetric shape turns 
out to bep*=(0.5M*)/(4/37!7-3) ~ 8.2xlO 9 M /fcpc- 3 . 

In both cases, given the estimations of the effective radius 
and the stellar mass of ~ 5xlO IO M , the source appears to be 
ultradense if compared with the local mass-size relation. In case 
(1) the stellar mass density should be considered an upper limit 
if the r e derived from the 1400A rest-frame is a fair estimate of 
the smallest size. Locally, on average, the r e is larger than 2 kpc 
(3 kpc) for early type (late ty pe) galaxies with the comparable 
stellar mass dShen et alj (120031) ) . 

Interestingly, the Lya modeling suggests (in all cases) a rela- 
tively high column density of the neutral gas, the Nhi turns out in 
the range 10 20,8-2 cm -2 . We further note that, assuming that the 
Schmidt law is valid at this redshift (Kennicutt 1998), adopting 
the observed area of 8.8 kpc 2 and two possible SFRs estimates 
(see Tabled, 30 and 100 M G yr~ l , the mass of the gas turns out to 
be 8 x 10 9 and 2 x 10 ll) M , respectively, which represents a sig- 
nificant fraction if compared to the stellar mass (~ 5xlO 1() M ). 

As noted by Buitrago et al. I (H008), massive ultradense 
spheroid observed at intermediate redshift ~ 1.5-3 and the glob- 
ular clusters have remarkably similar stellar densities (above 
lO lo M £pc~ 3 ), suggesting a similar origin. A massive ultradense 
galaxy at z~ 1.5-3 should form its stars very quickly in earlier 
epochs and in a high gas-density environment. In this sense the 
present source may represent a "precursor" of the ultradense 
spheroids recently discovered at redshift 1.5-3. 

5.3. Feedback in action ? 

The current burst of star formation may be caused by a previ- 
ous infall of gas and/or to a merger event (even though the UV 
morphology is quite regular). A vigorous wind of ~ 450 km/s 
is detected both from the observations (Lya and Niv] 1486 ve- 
locity offset) and from the Lya modeling (leaving all parame- 
ters free). As discussed above, the Lya emission arise from a 
compact region with an effective radius not larger than 0. 1 kpc 
(the PSF-deconvolved Re in the z'775 band is ~ 0.08 kpc physi- 
cal), a possible indication that the outflow of gas is in its initial 
phase of expansion from the central region. This particular phase 
of the galaxy evolution showing hot and massive stars and/or a 
low-luminosity AGN may be an efficient mechanism to blow the 
material out from the potential well of the galaxy, in some way 
influencing the subsequent star formation activity and/or the sur- 
rounding IGM. 

Wind propagation and escape is quite sensitive to the entrain- 
ment fraction and to the velocity of the wind itself. This occurs 
because the two primary forces limiting wind propagation are 
the galaxy's potential well and the ram pressure of the gas that 
must be swept up even if the wind is fast. Moreover, if entrain- 
ment is significant, then the mass over which the wind energy 
and momentum must be shared may be much greater. 



Therefore it is first useful to compare the escape velocity 
from the halo with the estimated w ind velocity. Following the 
calculation of iFerrara et al. I (l2000h . the escape velocity can be 
expressed as 

with p-\.65. The isothermal halo density profile is assumed 
(PH( r ) - Pc/[1 + ( r / r a) 2 ]), with an extension out to a radius 
? "200-'"/f = [3MH/47r(200p c ,,,)] 1 ^ 3 , defined as the radius within 
which the mean dark matter density is 200 times the critical den- 
sity p cr ii = 3Hq(1 + z) 3 /87rG at redshift z of the galaxy. The r# 
turns out to be -30 kpc at this redshift, assuming a halo mass of 
Mh - 1O 12 M . Under these assumptions, the escape velocity is 
v e =i 540kms~ l , about the same (or a bit larger) as the wind es- 
timated velocity from spectral features, ~450&m.!r I (or from the 
Lya modeling, ~500kms~ l ). It is possible that we are observ- 
ing the transport of material from the galaxy to the halo, which 
will remain confined. If we assume a slightly lower mass of the 
halo, e.g. M H - 5 x 1O U M (a factor 10 higher than the stellar 
mass), then the escape velocity turns out to be v e 3&0kms~ l 
and an r# =24 kpc. In this case, the velocity of the wind would 
be sufficient to escape the potential well of the halo. 

Indeed, from SPH simulations it appears that the main con- 
tributors to the metal enrichment of the low-density regions 
of the I GM are "smal l " galax ies with stellar masses below 
10 10 M o (Aguirre et al. (200M), 

and similar results have been 
obtained by other a uthors (e.g. [Oprjenheime & Dave I (2008), 
iBertone et al. I d2005l) ). In the present case the uncertainty on the 
halo mass prevents any clear conclusion. If we assume a value 
lower than 10 12 M Q , then the expanding material may reach char- 
acteristic distances (namely "stall radius") where the outflow 
ram press ure is balance by the I GM pressure up to few hundred 
kpes (e.g. Agui rre et aT71 ^OOlah t. 



6. Concluding remarks 

A peculiar galaxy belonging to the GOODS-S field has been 
discussed. The main observed features are the relatively strong 
nebular emission in the ultraviolet (Lya- and Niv] 1486) and 
the presence of the Balmer Break detected through the NIR 
VLT/ISAAC and Spitzer/IRAC data. Indeed, from the SED fit- 
ting with single and double stellar populations and the Lya mod- 
eling, it turns out that the source seems to have an evolved com- 
ponent with stellar mass of ~5xl0 10 M Q and age ~ 0.4 Gyrs, 
a young component with an age of ~ 0.01 Gyrs (contributing 
to ~ 1 % of the total stellar mass), and a star formation rate in 
the range of 30-200 M Q yr~ l . At present no evidence of com- 
mon "Niv] emitters" is observed in surveys of high redshift 
galaxies or quasars. However, there are rare cases in the liter- 
ature that show this line emission (together with other atomic 
transitions), spanning from a pure HII region source to a sub- 
class of low-luminosity quasars. In the first case, very hot and 
massive stars with low metallicity are required to produce the 
N iv] line; however, it is difficult to reproduce the signal mea- 
sured in the Spitzer/IRAC channels with a pure HII nebula, in 
particular at wavelengths beyond 4.5 pm, i.e. to reconcile the 
two observed facts: 1) the presence of a relatively evolved stellar 
population and 2) the low-metallicity environment needed if the 
N iv] emission arises from stellar photoionization. Alternatively, 
the low-luminosity quasar/AGN interpetation may explain the 
N iv] emission, the broad Lya component, and the properties of 
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the host galaxy discussed here, i.e., starforming, massive, and 
evolved galaxy. 

The limits on the size derived from the ACS/zgso an d VLT/Ks 
bands indicate that this object is denser than the local ones 
with similar mass, with a significant mass of the gas still in 
place (comparable to the stellar one). A relatively high nebular 
gas column density is also favored from the Lya line model- 
ing, Nffi>10 2l cm~ 2 , comparable to those found for the damped 
Lyman-alpha systems. The region emitting Lya photons is spa- 
tially compact, close to that at the continuum emission at 1400A, 
~ 0.1 kpc, in which a vigorous outflow (~ 450/500 km/s) has 
been measured from the spectrum and Lya modeling. The gas 
is expanding from this region, but given the uncertainty on the 
halo mass, it is dubious whether it will pollute the IGM to great 
distances. 

Such special objects are the key to understanding fundamen- 
tal passages in the formation and evolution of the galaxy pop- 
ulation. Future instruments will shed light on the nature of this 
interesting object, in particular, the JWST and the ELTs will give 
better and new constraints on the optical rest-frame morphology 
and nebular emission. 
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